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Abstract. Isospin violating hadronic decays of the rj and rj' mesons into Sir mesons are driven by a term 
in the QCD Lagrangian proportional to the mass difference of the d and u quarks. The source giving 
large yield of the mesons for such decay studies are pp interactions close to the respective kinematical 
thresholds. The most important physics background for 77, 77' — >■ tttttt is coming from direct three pion 
production reactions. In case of the rj meson the background for the decays is relatively low (« 10%). The 
purpose of this article is to provide an estimate of the direct pion production background for the rj' — > Stt 
decays. Using the inclusive data from COSY- 11 experiment we have extracted differential cross section for 
the pp — > pp-multipion production reactions with the invariant mass of the pions equal to the r;' meson 
mass and estimated an upper limit for the signal to background ratio for studies of the 77' — >■ tv'^tt'tv'^ 
decay. 

PACS. 13.60.Le Meson production - 13.75.-n Hadron-induced low- and intermediate-energy reactions 
and scattering (energy < 10 GeV) - 13.85.Lg Total cross sections - 25.40.-h Nucleon-induced reactions - 
29.20.Dh Storage rings 



1 Motivation 

1.1 Three pion decays of the rj and 77' mesons 

The r] and 77' decays into three pions violate isospin and 
occur only due to u and d quark mass difference. The decay 



width is sensitive to the mass difference: F, 



X 



ri(-q')- 



^0 • {fnd — mu)'^, where Iq term can be calculated in the 
isospin limit rud = ruu ■ The decays might provide a precise 
constraint for the light quark mass ratios [T]. 

In the case of the 77' meson the existence of the isospin 
conserving decays into three pseudoscalars (tttttj) implies 
that instead of the decay width for the isospin violating 
TTTTTT channel one can measure ratio: 



BR{r]' — > tttttt) 

BR{r]' -^ TTTTT]) ' 



(1) 



as it was proposed by Gross, Treiman and Wilczek [2|- 
Such measurement is self contained since it does not re- 
quire normalization for the partial decay width from other 
experiments. Additionally, simultaneous measurement of 
the two decays modes, with similar final states, ensures 
that many systematic uncertainties will cancel. The 77' — >■ 
TTTTTT decays provides also a very sensitive test of the Chiral 
Perturbation Theory (CHPT) framework [31|3] extensions 



to the 77' meson. Due to the large mass of the 77' meson, the 
decays are strongly influenced by light vector and scalar 
meson resonances. Those decays cannot be studied using 
standard CHPT methods. An elegant method for account- 
ing for the nonperturbative effects in two pseudoscalar 
meson interactions is given by unitarization procedure of 
the one loop CHPT result. For example in the theoreti- 
cal studies of the 77 and 77' meson decays the rescattering 
of any pair of the pseudoscalar mesons is described by 
Bethe-Salpether equations "W^^. The parameters of the 
interactions are obtained by fits to the pseudoscalar scat- 
tering data. Predictions for many 77 and 77' decays were 
given using the above technique within so called chiral 
unitary approach [8]. 



Experimentally determined value of the branching ra- 
tio of the 77' -^ 3tt° decay is (1.56±0.26)xl0-3 [9]. The 
decay into Tr"'"7T~TT° was observed in 2009 for the first time 
by the CLEO collaboration and the branching ratio was 
determined to (37t" ± 4) x 10""* 10]. The result is in 
strong disagreement with the value 10~^ predicted within 
the chiral unitary approach. Much more improvement on 
the theory and experiment side is needed to understand 
the three pion decays of the 77' meson. 
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Fig. 1. Compilation of the data for the pp — )■ pprj' reaction 
cross section from COSY- 11, DISTO and SPESIII measure- 
ments [22ll23U24U26ll27| . The solid line is parameterization of 
the data using formula ^. 



1.2 Experiments using pp — >■ ppr]',ri reaction 

The mesons for the decay studies are produced in jp [TTJ 
[in, pp Bun, n-p [UlIIS], pd ^MMM or e+e- ^U^^ 
interactions. For the studies at hght ion storage rings as 
COSY the pp — ^ pprj' reaction close to threshold seems to 
be most promising. The cross section for pp — > ppr]' reac- 
tion was measured by the COSY-11 [111[2S1[M1[^ , SPE- 
SIII [26] and DISTO [27l coUaborations. In Fig. [U the ex- 
perimental data are compared to the analytical parame- 
terization derived by Faldt and Wilkin [28 , 25] which takes 
into account final state interaction of the protons: 



'(Q) = c- 



Q' 



1 



rUpPLAB I 



1 



where Q denotes the excess energy, plab beam momen- 
tum and rup proton mass. In comparison to the pp inter- 
action the p — r]' interaction is negligible [301 ■ The C and 
e free parameters were determined by a fit to the experi- 
mental data [S]: e = 0.62 ±0.13 MeV and C = 42±7 mb. 

Experience from studies of the rj — > tttttt decays with 
rj produced in pp — > pprj reaction at beam energies 1.30- 
1.45 GeV carried out by the CELSIUS/WASA collabora- 
tion shows that background from direct three pion pro- 
duction is about 15% for the tt+tt^tt" and about 5% for 
the 37r° channel [5^I33| . This allows for precise study of 
the 77 decays providing a large number of events is col- 
lected. The production cross section for the r] mesons in 
pp colisions [Mll551l5Bll37l[T3] is about 30 times larger than 
the cross section for the 77' meson [^123,24,25,26,27 at 
the corresponding excess energies. At the same time the 
total cross section for the direct three pion production in- 
creases about two orders of magnitude between 77 and 77' 
production thresholds (Fig. [5]). For the pp — >■ pp7r+7r~7r" 
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Fig. 2. Total cross section for three pion production: data and 
parameterizations. The data are from '43',' 44, 45, 46, 32'. The 
parameterizations are from Bystricky et al. 47 and Faldt and 
Wilkin [28II29| . The kinetic beam energy at thresholds for the 
rj and rj' production is equal to 1.255 GeV and 2.404 GeV 
respectively. 



reaction total cross section there are only three exper- 
imental points in the beam kinetic energy range up to 
3 GeV. The cross section for the pp — )■ pp7r°7r°7r° reac- 
tion near the rj' meson production threshold was not mea- 
sured at all. However, based on the statistical model [551 
l5^Hn] and the isobar model [CTH^ one expects the the 
pp — >■ ppTT^TT^Tr" cross section to be 6— 10 times lower than 
for the pp — >■ ppn'^ TT~ TT^ reaction. This is in agreement 
with an extrapolation of the CELSIUS/WASA measure- 
ment of the both reactions close to 77 meson production 
threshold. For the estimate of the background for the three 
pion decays of the rj' mesons instead of the total cross sec- 
tion the relevant quantity is the value of the differential 
cross section for the invariant masses of the pions in the 
range of the 77' meson mass. This quantity was not mea- 
sured and in the present paper we will provide an estimate 
for the upper limit of the background using the COSY-11 
data where only outgoing protons were registered. 



2 General considerations 

Let us consider an example analysis chain leading to a 
selection of the 77' — >■ 7r'''7r~7r° decay from the pp — > ppr]' 
reaction. In the first step all tracks are reconstructed and 
particles identified. The events containing two protons 
from the production process, two charged pions and two 
photons are selected. Now one can apply energy-momentum 
conservation and select only events consistent with pp — ^ 
pprj' — )• ppTT^TT^n'^ — >■ pp7r"'"7r~77 reaction hypothesis. 
This procedure can be most generally implemented by 
kinematic fitting but one can use also some other method. 
In the end the selection could be represented by a region 
in some control variable p.. For example p could be miss- 
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ing mass squared or x^ value of the kinematic fit. Witlrin 
the selected region, in addition to the signal, a contami- 
nation from the background events originating from reac- 
tions which have similar final states is unavoidable. The 
most important physics background channels for the dis- 
cussed case is the direct pp — ^ pptt^tt'tt^ reaction and 
other rj' decays like rj' — >■ tt'^'tt"?] — >■ 7r+7r~77. Candidates 
for the identification variable (in addition to x^ of the 
kinematic fit) are in these cases respectively: the missing 
mass of the two protons and the invariant mass of the 
two photons. In this article we focus on the first case: 
estimate the signal to background ratio and its implica- 
tions for the statistical uncertainty of the extraction of 
BR{ri' -^ 7r+7r"7r°) value. 

The signal to background ratio, Ns/Nb, can be writ- 
ten as: 



Ns_ 
Nb " 
where the factors are: 



arji ■ BR ■ £s ■ C 
A^,- pb ■ eb ■ C 



(3) 



1. Gri' - the total cross section for the production reaction 
(here for pp — > pprj'), 

2. pb - the differential cross section for the direct tt+tt^tt^ 
production with the pions invariant mass equal to the 
mass of the rj' meson: 



duf 



PB 



d/i 



(4) 



6. 



£St £b ~ acceptances and reconstruction efficiencies for 
the signal and the background, 

Ap - a range of the missing mass used for the extrac- 
tion of the signal it depends on the detector resolution, 
BR - the measured branching ratio of the rj' -^ tt^tt~tt^ 
decay, 
C - stands for the integrated luminosity. 

The Ns/Nb ratio depends on the beam energy through 
o't]', Pb, the missing mass resolution and the detection 
efficiencies. Hereafter we derive the energy dependence of 
these quantities. 



3 Background estimate 

The value of the pB cross section should be determined 
from the 7r+7r~7r° invariant mass distributions of the di- 
rect pion production reaction. However, there is no data 
at beam energies near the 77' meson threshold. Therefore, 
we estimate an upper limit for the pB by re-evaluating the 
available missing mass spectra of the pp — > ppX reaction 
determined by the COSY-11 collaboration at several beam 
energies near the the 77' threshold [22,23,24,25 . In Fig. [3] 
an example of COSY-11 reconstructed missing mass dis- 
tribution at the pp — ?> ppr]' excess energy, Q, of 15.5 MeV 
is shown [511 . 

Let us consider a measurement of the two protons from 
the pp — >■ ppX reaction, where there are no constrains on 
other outgoing particles. In a first approximation the ratio 
of the number of the background events in a slice under 
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Fig. 3. An example of the missing mass distribution for the 
pp — !■ ppX reaction from the COSY-11 measurements at Q = 
15.5 MeV ^I]. S denotes the signal originating from the r;' 
meson production and B indicates the background under the 
peak. The continuum originates from direct production of two, 
three and more mesons. It provides a conservative upper limit 
for the background from the pp -^ pp7r"'"7r~7r° reaction. 



the 7]' peak in the missing mass spectrum to the number of 
events in the peak does not depend on the detector accep- 
tance for the protons. This assumption is valid for example 
if the rj' meson and the multipion reactions are simulated 
according to phase space or even if an universal final state 
interaction between protons is introduced [48ll49] . The as- 
sumption may break for example if there will be a signif- 
icant difference in four momentum transfer distributions 
between rj' and multimeson production. However, the de- 
pendence of the production amplitude on the momentum 
transfer is very weak near the kinematical threshold. 

The differential cross section pB for the background 
originating from all multimeson channels was determined 
from the COSY-11 data according to the formula: 



Pb{Q) = 



NB{Q)<7rr{Q) 
NsiQ) Ap ■■ 



(5) 



which is derived from the following expressions for the 
iVs(Q) and iVB(Q): 

Ns{Q)=cT^'{Q)-e{7nr,,,Q)-C, (6) 

NB{Q)=PB{Q)-Ap-e{mn,,Q)-C, (7) 

where Ng stands for the number of the observed events in 
the rj' peak, Ns number of the background events in the 
Ap slice under the ry' signal, ct^/ denotes the total cross 
section for the pp — >■ pprj' reaction described according to 
analitical formula from Eql^l e denotes combined accep- 
tance and detection efficiency of the COSY-11 detector, 
which, in a very good approximation, depends only on 
the mass of the produced system and on the excess en- 
ergy P51l49) ■ C indicates the integrated luminosity. The 
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Table 1. Differential cross section pB for multimeson produc- 
tion in proton-proton collisions as a function of the pp -^ pprj' 
reaction excess energy Q. The pB values were extracted from 
the experimental data |22II23II24|[^ using Eq. [S] 
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Fig. 4. Inclusive differential cross section for multimeson pro- 
duction derived from the COSY- 11 data [22ll23U24ll25] . Statisti- 
cal and systematic errors are seperated by the horizontal bars. 
The superimposed line shows the function of Eqn. |8] fitted to 
the COSY- 11 data. 



width of the Afi slice was selected to contain nearly 100% 
events of the signal peak. The derived values of pB as a 
function of Q with corresponding statistical and system- 
atic errors are given in Table [T] and are shown in Fig. Ul 
The systematic uncertainties are discussed in [22l[23l[24l 
[25]. 

The pb dependence on the excitation energy is well 
described by the following parameterization: 



PB = aiQ/QoY 



(8) 



where Qo — 1 MeV is the normalization factor and a 
and (3 are the parameters. The fit gives a = 0.64 ± 0.14 
nb/MeV and (3 = 1.662 ± 0.081. 



4 Statistical uncertainty of the BR 
measurement 

Assuming that the shape of the background is known and 
is not correlated with the signal the relative statistical 
uncertainty of the branching ratio can be expressed as: 



ajBR) _ a{Ns) _ y/Ns + Nb 
BR " Ns ^ Ns 

In our case Ns <C Nb and we have: 

a{BR) _ y/NB 
BR ^ Ns 



(9) 



(10) 



Taking into account Eq.|3]and Eg.lTUlone gets the formula 
for the statistical uncertainty: 



a{BR) < 



yj Pb ■ Ap.- eb 1 



£s 



£ 



(11) 



From the derived expression one sees that improved tag- 
ging resolution {Ap decreased) helps only if the detection 
efficiency is not worsened. 

The integrated luminosity H can be determined from 
the simultaneously measured decay 77' — >■ 7r+7r~77 — >■ 7r"'"7r~77 
with well established branching ratio. Therefore the sta- 
tistical uncertainty of the luminosity determination can be 
neglected and many contributions to the systematic un- 
certainty in the ratio of signal to the monitoring events 
will cancel. 



5 Feasibility for a large acceptance detector 

As an example application of the extracted pB value and 
the formulas derived in the previous sections we will con- 
sider a determination of the branching ratio for the 77' — ^ 
7r"'"7r~7r° decay using a large acceptance detector. The 
large acceptance is necessary for efficient identification of 



all outgoing particles. After selection of the pp7r+7r~7r ' fi- 
nal state the direct tree pion production and the 77' — > 
TT+Tr^Tr*^ decay can be distinguished with the best preci- 
sion using the missing mass of the two forward emitted 
protons. 

For the calculations of the missing mass resolution 
a typical beam and target parameters available at the 
Cooler Synchrotron COSY [50] were assumed: beam mo- 
mentum spread Ap/p « 10^"^ (FWHM), perpendicular 
beam profiles: horizontal crx=2 mm, vertical aY=5 mm 
[5T] . As a target a hydrogen stream in a cylinder with di- 
ameter of 2.5 mm was used. Effective energy resolution of 
the forward scattered protons from the reaction pp — >■ ^^77' 
registered using plastic scintillators is typically in the or- 
der of few percent (3%). 

The determined Q dependence oia{BR)/BR is shown 
in Fig. [HI assuming one week experiment with luminosity 
of L = lO^^cm-^s-i and value of BR{q' -^ tt+tt-ttO) = 
0.37% [TU]. The optimum is reached for the excess ener- 
gies between 50 and 100 MeV. This is a general conclusion 
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Fig. 5. The relative accuracy of the upper hmit of the 
BR{rj' —J- 7r"'"7r~7r'') as a function of the excess energy Q for 
the pp — >■ pprj' reaction. 



Fig. 7. The signal to background ratio S:B calculated taking 
into account the natural width of the rj' meson. 
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Fig. 6. The relative accuracy of the upper limit of the branch- 
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for a large acceptance detection systems. The statistical 
uncertainty of the branching ratio improves with time of 
the measurement as l/yt. The dependence for the beam 
momentum of pbeam — 3.45 GeV/c corresponding to the 
excess energy Q = 75 MeV is shown in Fig. [51 The plot in- 
dicates that for the BR equal 0.37% to achieve the relative 
accuracy of 10% would require at least two month exper- 
iment. However signal to backround ratio at this energy 
will be only about 10~^ what puts extreme requirements 
for the understanding of the systematic effects. Therefore 
the other strategy for the experiment would be to find a 



compromise between the statistical and systematic uncer- 
tanities by going to lower excitation energies where the 
signal to background ratio increases (FigjT]). 

An additional source of background, not discussed here, 
comes from other decays of rj' involving similar particles: 
rj' — >■ Tv^Tv^T] and rj' — > ujj. This background cannot be 
suppressed using the missing mass method and the invari- 
ant masses of the decay products should be used instead. 



6 Summary 

Using the COSY-11 data for the pp — )■ ppX reaction near 
the rj' meson production threshold we extracted an up- 
per limit for the background from the pp — > pptv^tt^tt^ 
reaction for rf — )• 7r+7r~7r*' decay studies. This and the 
parameterization of the total cross section energy depen- 
dence for the pp — >■ ppr]' reaction permit us to estimate 
that two months of the beam time with average luminos- 
ity of lO^^cm^^s^^ would be sufficient to reach statistical 
accuracy of the previous experiments for the studies of 
77' — >■ 7r+7r~7r° using a large acceptance detector. However 
since signal to background ratio for the optimal energy is 
only about 10"'^ one may expect much larger systematic 
uncertainty. Since even taking the natural width of the rj' 
meson the S : B ratio is about 10~^, this points to inherent 
limitations of the pp — >■ pprj' reaction for the rj' — ^ 7r+7r~7r'' 
decay shown in Fig. [T] Situation is expected to be at least 
one order of magnitude better for the rj' — )■ tt^tv^tv^ case. 
Also "not rare" decays 77' 
pp — >■ ppf]' reaction. 



TTTTTy can be studied in the 
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